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Abstract 



We report grand canonical Monte Carlo simulations of the critical point 
properties of homopolymers within the Bond Fluctuation model. By employ- 
ing Configurational Bias Monte Carlo methods, chain lengths of up to N = 60 
monomers could be studied. For each chain length investigated, the critical 
point parameters were determined by matching the ordering operator distri- 
bution function to its universal fixed-point Ising form. Histogram reweighting 
methods were employed to increase the efficiency of this procedure. The re- 
sults indicate that the scaling of the critical temperature with chain length is 
relatively well described by Flory theory, i.e. G — T c ~ iV~ ' 5 . The critical 
volume fraction, on the other hand, was found to scale like <f> c ~ iV~ - 37 , in 
clear disagreement with the Flory theory prediction (f) c ~ iV~ ' 5 , but in good 
agreement with experiment. Measurements of the chain length dependence 
of the end-to-end distance indicate that the chains are not collapsed at the 
critical point. 



PACS numbers 61.25.Hq, 64.70.Fx, 05.70.Jk 



I. INTRODUCTION AND OVERVIEW 



When long flexible polymers are dissolved in a bad solvent there exists a critical tem- 
perature T C (N) of unmixing slightly beneath the B-temperature (figure [1]). At this critical 
temperature, the system phase separates into a very dilute (solvent rich) solution of collapsed 
chains and a semidilute (polymer rich) solution. The process is qualitatively described by 
the mean field theory of Flory which predicts simple power laws for the chain length 
(A) dependences of T C (N) and the corresponding critical volume fraction (f) c (N): 

T C (N) = Q/(1 + 1/VN) 2 (1) 
^0-26/^, A^oo 



C (A) = 1/(1 + VN) (2) 
ml/VN, A^oo. 

Another power law is predicted for the shape of the coexistence curve near T C (N): 

€L - *2L = 2B(N)t^ t=1- T/T c (A), (3) 

with a critical order parameter exponent j3 and a chain length dependent critical amplitude 
B(N) given by 

13 = f3 MF = 1/2, B(N) oc A- 1 / 4 . (4) 

Further power laws describe the intensity of critical scattering, the associated correlation 
lengths and the interfacial tension etc but will not be considered here. 

Notwithstanding the qualitative correctness of the Flory theory in predicting a phase 
separation, it should be emphasised that the exponent (5 = /3mf = 1/2 in equation f|, as well 
as the powers of A in equations 0-§] are mean field results, and thus cannot be expected 
to be quantitatively correct. More generally one expects that (we follow the notation of a 
recent experimental study ||) 

5(A)ocA- xl , (5) 
C (A) oc A~ X2 , 
6-T c (A) oc N~ X3 , 

where the mean field values of the exponents defined in equation [| are 

*r = 1/4, (6) 

It is an interesting open question to ask what are the correct values of these exponents. 
While it is generally accepted from the "universality principle" 0, as well as experimental 
findings P,|7|- |T3[1 , that the phase separation of polymer solutions falls in the same universality 
class as the three dimensional Ising model, so that || 



1 



(3 « 0.325, 



(7) 



the theoretical understanding of the exponents X\,x 2 ,X3 in equation |5| is rather limited. 
Experimental data have yielded the estimates P,|7j-|l4"l 

x t « 0.23 - 0.34, (8) 

x 2 w 0.38 

x 3 w 0.47- 0.50. 



However, theoretical estimates for these exponents are still controversial. De Gennes fl5 
suggested that in the limit of large N, one has the same scaling behaviour as in mean field 
theory, i.e. the coexistence curve scales as 



-6® = — 

coex "coex 



, n <KvNt), (9) 

Since <p(z) must behave for small argument z as <j>(z) oc z 13 , this yields 

Xl = (1 - p)/2 » 0.34, (10) 



which is roughly compatible with experiments. However, the scaling with v N in equation |9] 
implies that x 2 = 1/2, which clearly disagrees with equation |8|. Muthukumar [IB] on the 



other hand, suggested that in a limit where ternary interactions are important, one should 
have different exponents, namely 

Xl =x 3 = 2/9, x 2 = l/3. (11) 

Subsequently this problem has received further attention in the literature [§,|17]-|20| . Re- 



call that the scaling structure in equation |D| can be justified in terms of a Ginzburg criterion 
21j| if one assumes that the chain linear dimensions are ideal ]15| . It then follows (remember- 



ing that the chain gyration radius enters as a critical amplitude prefactor in the mean-field 
power laws of the correlation length of the monomer density fluctuations), that the critical 
density c (iV) coincides (up to a universal prefactor) with the onset of the "semi-dilute" 
regime, where chains overlap significantly ||15|| . This assumption is plausible because of the 
vicinity to the 9-state (T = 6, <fi — > 0), where chains indeed behave ideally and the gyration 
radius scales as 

i? 5 ociV 1/2 , T = 6, 0->O, N^oo. (12) 
However the fact that for T < and (f) — > 0, N — > oo chains are collapsed: 

i? 5 ociV 1/3 , T<9, 0^0, iV^oo, (13) 
implies that one does not really know how R g scales with iV at the critical point. Therefore 



it is tempting to generalize the scaling ansatz R] as follows [19 



W -rhW =J_0(JV«*r). (14) 
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Equation H is, of course, still consistent with the behaviour of the coexistence curve at 



fixed r in the limit N — > oo [[- 
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coex 



coex 



-(l 

2 V 
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(15) 



if 



oo 



\z. Since for small z = N X4 t, the scaling function must behave as <j>(z) oc z- 



P 



in order to comply with equations |3] and [7], we conclude that 
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id) 

coex 



^4(1-/3) a . i =^(1-/3). 



(16) 



From renormalization group arguments, Cherayil |L9| has suggested that the exponents x 2 
and £3 can be expressed in terms of the new exponent x± as 



x 2 



X4, X3 — X4.. 



(17) 



This theory, however, does not yield a prediction for X4 itself, and to fit some experimen- 
tal data it was assumed that £4 = 0.62, x 2 = 0.38 Kholodenko and Qian |Tj| have 
presented arguments that the exponent x 2 is not even a universal quantity. If the scaling 

hold, this would imply that x±, x 2 and x% are all system 



relations of Cherayil (equation [TBlJTi 
specific quantities, depending upon the material under consideration! Finally, we note that 
Muthukumar's result, equation 11, disagrees with the above scaling relation X\ = 23(1 — (3), 



and thus the theoretical situation is clearly somewhat confusing. 

In view of these problems and the difficulties of extracting all relevant information from 
experiments (one not only wishes to check the relations of equation ^ but also seeks to 
clarify how the chain span scales with N at criticality), study of this problem by Monte- 
Carlo computer simulations techniques |Z2 is highly desirable. In fact there has been some 
previous work on this problem which considered the vapour-liquid phase diagram of alkane 
chains [23| and coarse-grained off-lattice polymer models (see e.g. [p4] ,^ 
work of reference [p3| 



However the 



[see e.g. 

considers the problem of estimating absolute values of T C (N) and 
C (N) for a chemically realistic model of alkanes for small N and does not address the 
universal properties of the limit iV — ► 00. The Gibbs ensemble Monte-Carlo method of 



Panagiotopoulos |26| , [27|] allows an efficient estimation of the coexistence curve well below 
the critical point, but a precise estimation of critical point parameters is difficult in this 
framework. 

An alternative approach for estimating critical point properties from simulations is based 
on finite-size scaling f28| , |2"9"|| . This approach has been very successful for both symmetrical 
31 1 and asymmetrical |32| polymer mixtures in conjunction with the bond-fluctuation lattice 



model [35] and semi-grand canonical ensemble simulation techniques p0 |. These studies also 



relied on the use of histogram reweighting |33] and (in the asymmetric case) recent advances 
in disentangling order parameter and energy fluctuations near criticality in a finite-size 



scaling context |34| . 

In the present work we attempt to apply a related approach to study the liquid- vapour 
critical point of homopolymers within the Bond Fluctuation model. This problem is, how- 
ever, somewhat more intricate than that of polymer mixtures since one must employ the 
grand canonical ensemble (GCE) [|34| in order to effectively deal with the strong near-critical 
density fluctuations. As is well known, GCE simulations for chain molecules are extremely 
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difficult, since the insertion probability for a polymer chain into a many chain systems is 
vanishingly small [^2,[36-38|]. For chains that are not too long (and/or systems that are not 
too dense), this problem can be eased by the Configurational Bias Monte-Carlo (CBMC) 
Method |36"H38|1. In the present paper we combine CBMC with histogram reweighting and 
a finite-size scaling analysis in its form extended to asymmetric systems [^,0. By this 
special combination of recent techniques (which will be briefly reviewed in section ||) we are 
able, for the first time to obtain accurate results, both for 4>(N) and T C (N) up to N = 60 
effective monomers. Since the effective bond in the bond-fluctuation model can be thought 
of as corresponding to 3 to 5 chemical bonds (when a mapping to chemically realistic chain 
molecules is attempted [0), our simulations thus correspond to a degree of polymerization 
up to a few hundred chemical bonds along the chain backbone. 

Section [III] then presents our results, including a highly precise estimation of the 
temperature from an analysis of the gyration radius of simple isolated chains. We obtain 
both the location of the critical point in the (T, 0) plane as a function of chain length and, 
for the first time, the associated dependence of the chain span. In section [IV| we discuss 
our results and compare them to the theoretical ideas sketched above. We obtain very 
good agreement with experiment, but as in the latter the need to study much longer chains 
is clearly apparent to definitively clarify the true asymptotic behaviour for chain lengths 
N -> oo. 



II. ALGORITHMIC AND COMPUTATIONAL ASPECTS 

The bond-fluctuation model (BFM) studied in this paper is a coarse-grained lattice-based 
polymer model that combines computational tractability with the important qualitative fea- 
tures of real polymers, namely monomer excluded volume, monomer connectivity and short 
range interactions. Within the framework of the model, each monomer occupies a whole 
unit cell of a 3D periodic simple cubic lattice. Neighboring monomers along the polymer 
chains are connected via one of 108 possible bond vectors. These bond vectors provide for 
a total of 5 different bond lengths and 87 different bond angles. Thermal interactions are 
catered for by a short range inter-monomer potential. The cutoff range of this potential was 
set at r m = v^6 (in units of the lattice spacing), a choice which ensures that the first peak 
of the correlation function is encompassed by the range of the potential. We note also, that 
within our model, solvent molecules are not modelled explicitely, rather their role is play by 
vacant lattice sites. Further details concerning the BFM can be found in reference . 

To implement a grand canonical ensemble simulation of the BFM, the Configurational 
Bias Monte Carlo (CBMC) method was employed |3"E|-|3"8"|. The CBMC scheme utilizes a 



biased insertion method to 'grow' a polymer into the system in a stepwise fashion, each 
successive step being chosen so as to avoid excluded volume where possible. For brevity we 
shall merely outline the GCE implementation of this CBMC method and refer the reader 
to reference for a fuller description. 

Within the GCE scheme there are two complementary types of moves, insertion attempts 
and deletion attempts, both of which are made with equal frequency. An insertion move 
first involves attempting to grow a candidate polymer into the system. The basic strategy 
for achieving this is to insert successive monomers of the chain into the system one by one. 
The position of each successive monomer is chosen probabilistically from the set of 108 
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possible BFM bond vectors emanating from the previously inserted monomer. The selection 
probability for each of the possible monomer positions is weighted by its Boltzmann factor, 
effectively biasing the choice in favour of low energy chain configurations. In order to keep 
track of the accumulated bias, a book keeping scheme is maintained. Once a candidate chain 
has been successfully grown, it is submitted to a Monte-Carlo lottery to decide whether or not 
it is to be accepted. The total chain construction bias is compensated for in the acceptance 
probability, thereby ensuring that detailed balance is obeyed. 

For chain deletion moves, one chooses a chain at random from those in the system and 
'reconstructs' its bias by examining the alternative growth scenarios at each step of the 
chain. The candidate chain for deletion is also submitted to a Monte-Carlo lottery to decide 
whether the proposed deletion should take place. As with the insertion lottery, the chain bias 
is taken into account in the deletion probability. The chemical potential, fj,, which controls 
the system chain density, also enters into the acceptance probability for both insertion and 
deletion. 

The principal observables measured in the course of the simulations were the monomeric 
volume fraction: 

cf) = 8nN/V (18) 

and the dimensionless energy density: 

u = %w- l $({r})/V (19) 

where n is the number of chains, $({r}) is the configurational energy, w is the depth of 
the square well interaction potential (so that T = w _1 ) and V is the system volume. Here 
the factor of 8 derives from the number of lattice sites occupied by one monomer in the 
BFM. Measurements of <fi and u were performed at intervals of 500 chain insertion attempts 
and accumulated in the joint histogram pl{4>,u). The final histograms comprised some 10 5 
entries. Also measured were the distributions of the chain radius of gyration and the chain 
end-to-end distances. 

Using the GCE algorithm, chains of lengths N = 10, 20, 40, 60 were studied. For the 
N = 10 and N = 20 system size V = 40 3 and V = 50 3 were employed, while for N = 40 
and N = 60 chain lengths only the V = 50 3 was studied. Unfortunately it was not possible 
to study chains longer than N = 60 since the acceptance rate for chain insertions falls 
exponentially with increasing N and volume fraction cf). This problem is illustrated in 
figure & where we plot the acceptance rate for a number of chain lengths as a function of 
the monomeric volume fraction. One sees for example, that for N = 80 the acceptance 
rate is too low to provide reliable statistics within reasonable run times. Indeed, even for 
our longest chain length N = 60, extremely long runs were required to gather adequate 
statistics. 

Having outlined our model and simulation technique we now turn to a brief description 
of our data analysis methods. As mentioned in the introduction, finite-size scaling (FSS) 
methods are an indispensable tool for the proper treatment of critical behaviour, serving as 
they do to provide estimates infinite-volume critical properties from simulations of finite- 
size systems. The FSS methods we shall employ here are especially tailored to fluid systems 
and have been described in detail elsewhere |34j. The basic idea is to exploit the Ising 
character of the polymer liquid- vapor critical point to accurately locate the critical point. 



5 



This is done by observing that precisely at criticality the distributions of certain readily 
measurable observables assume scale-invariant universal forms. The particular universal 
scaling form on which we shall focus, is the distribution of the ordering scaling operator 
Pl(M) . For the Ising model, the special symmetry between the coexisting phases implies 
M. —>■ m (the magnetisation). The critical point form of Pl(^) is independently known 
from extensive studies of large Ising lattices p9[ . For fluids, however, the lack of symmetry 
between the coexisting phases implies Pl|| that the ordering operator is a linear combination 
of the fluid density and energy density i.e. M. — > <p + su i where s is a system specific 'field 
mixing' parameter that controls the strength of the coupling between the density and energy 
fluctuations. 

Thus, in principle, one is able to accurately locate the critical point of a fluid system 
simply by tuning the T, /i and s until the measured form of Pl{M) matches the known fixed 
point Ising form. In performing this task, the histogram reweighting method |33| , |3l| can be 
of great assistance. This technique allows one to generate estimated histograms pi((f),u) for 
values of the control parameters T and \i other than those at which the simulations were 
actually performed. Such extrapolations are generally very reliable in the neighborhood of 
the critical point, due to the large critical fluctuations ||33|| . In what follows we shall detail 
the application of all these techniques to the problem of determining the liquid- vapor critical 
point parameters of our polymer model. 



III. PROCEDURE AND RESULTS 

The first task undertaken was a determination of the 0-temperature for our model, 
knowledge of which is a prerequisite for studying the scaling of — T C (N) (cf. equation ||). 
To achieve this the gyration radius R g of single chains was studied as a function of tempera- 
ture and chain length. From equations one sees that precisely at the 0-temperature, and 
modulo corrections to scaling, R 2 g /N should be independent of N. Extensive simulations 
were therefore carried out for single chains of length N = 64, 80, 100 and 150 at a tempera- 
ture T = 2.0. The full temperature dependence of R g for each chain length was subsequently 
obtained by histogram reweighting. This involves recording the joint histogram of the gy- 
ration radius and conformational energy of each configuration generated. The histogram 
for other temperatures may then be obtained by reweighting the Boltzmann factor for each 
histogram entry in the manner described in reference [33|. Figure ^| show the resulting curves 



of Rg/N(T), which exhibit a very precise intersection point at T = 2.020(3), a value that 
we therefore adopt as our estimate of the 0-temperature. 

In general for fluid systems, the coexistence curve is not known a-priori and must there- 
fore be identified empirically as a prelude to locating the critical point itself. In the following 
we exemplify the general strategy for determining the critical parameters by considering the 
case of the N = 20 system. 

Initially a temperature of T = 1.75 somewhat beneath the temperature was chosen 
and the approximate value of the coexistence chemical potential was determined by tuning 
/i until Pl{4>) exhibited a double peaked structure. A long run was then carried out at this 
near-coexistence /i value, in which the histogram of Pl{4>, u) was accumulated. A histogram 
extrapolation based on this data was then used to extrapolate along the coexistence curve 
using the equal peak-weight criterion for Pl{4>) In this way a sizeable portion of the 
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the near-critical coexistence curve (and its analytic extension ||34|| ) could be located. Repre- 
sentative forms of the density distributions along this line of pseudo-coexistence are shown 
in figure ffl(a), while the positions of the density peaks are shown in figure f|(b). 

To locate the critical point along the line of phase coexistence, we utilized the universal 
matching condition for the operator distributions Pl{M) ■ Again applying the histogram 
reweighting technique, the temperature, chemical potential and field mixing parameter s 
were tuned until the form of Pl{M) most accurately matched the universal fixed point Ising 
form fi^^x). The results of performing this procedure are shown in figure [5] for the V = 40 3 
and V = 50 3 system sizes. Given that these systems contain on average less than 100 polymer 
chains at criticality, the quality of the data collapse is remarkable. The mappings shown 
were effected for a choice of the parameters T c = 1.791(5), ji c = —5.16425(2), s = —0.135(4) 
where we have defined /i to be the chemical potential per monomer. The associated estimates 
for the critical volume fraction is <p c = 0.198(5). We note that this estimate of T c is rather 
less than what would be obtained, were one simply to extrapolate the liquid and vapour 
densities to the point at which they merge [figure £|(b)]. Thus our results further emphasise 
the finite-size errors that can arise using such a procedure. 

We have also attempted to fit the coexistence curve data of figure (|(b), using fits of the 
Ising form: fa = C +O.745(T C -T)+O.55(T C -T) a324 , <\> v = C +O.745(T C -T)-O.55(T C -T) a324 , 
with T c and C assigned the values obtained above. For the Lennard- Jones fluid |4*2|| , a fit 
of this type gave a good description of the coexistence data for temperatures down to 0.8T C . 
Here, one sees that a reasonable fit can be obtained, but only for temperature within some 
3% of T c . This would seem to suggest that the asymptotic Ising region in the polymer system 
is much smaller than in the Lennard- Jones system. 

With regard to our estimates for the critical parameters, it should be emphasized that 
they are subject to errors arising both from corrections to scaling, as well as field mixing 
effects in the case of the critical volume fraction [34] . While in principal one can also correct 
for these effects, if one has access to a sufficient range of system sizes and ample statistics 
(see eg. ^2|), the computational difficulties of the present problem preclude such an analysis. 
Indeed, for the N = 60 system studied here, only one system size, V = 50 was employed, this 
being the largest that could reasonably be tackled. Smaller system sizes were not studied 
since these would contain so few chains at criticality as to be excessively influenced by 
corrections to scaling. The matching to the universal form for the largest size available was 
then the only guide to the location of the critical point. Notwithstanding these problems 
however, we feel on the basis of our experience of corrections to scaling in other systems |4"2" 



that the quoted uncertainties generously encompass the infinite-volume critical parameters. 

The procedure described above was repeated for the other chain lengths studied, allowing 
estimates to be obtained for T C (N) and (f> c (N). However the computational difficulty became 
progressively greater as N increased, making the accumulation of good statistics problematic. 
As a result it was not possible to perform a reliable histogram extrapolation away from 
the critical point into the sub-critical two-phase region and thus no information on the N 
dependence of the critical amplitude prefactor featuring in equation |5| could be obtained. An 
additional hindrance to probing the subcritical coexistence region is that the BFM appears 
to be unable to support a liquid phase for volume fractions ^ 0.6, instead collapsing into 
an amorphous crystal. This artifact is traceable to the limited conformational entropy of 
our lattice-based chains, and has also been observed in a previous study of tethered chains 
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using the same model f43| . 

The results for T C (N) and C (-^O are plotted against iV -1 in figures || and ^respectively. 
For T C (N) we find that the data can be well fitted by a Flory-type formula of the form 
T c = O + c^iV -0 ' 5 + c^iV -1 (where the iV _1 term can be thought of as a free-end correction 
and where we assume that UmN-*ooT c = &), although fits of the form T C (N) = O + aiN~ X3 
yield a comparable fit quality for values of 23 in the range x% = 0.46-0.53. For <f)(N) we 
have performed a fit of the form c (iV) = (61 + b2N X2 )~ 1 , and obtain xi = 0.37(2). 

Finally we have considered the N dependence of the average squared end-to-end distance, 
R 2 ee (N) at criticality. This quantity has been conjectured to scale as 



RUN) oc N» (20) 

It was observed that the scaling of C with N could be explained if v' = 0.46, and we would 
like to check this conjecture. Our results are plotted in figure |8|. Despite the very limited 
number of data points, we have attempted to fit this data to the form Eq. ^(J, with the result 
R? ee (N) oc jV L11 ( 4 ). This finding that v' > 0.5 is also supported by a study of the distribution 
function of end-to-end distances p(r) at the estimated critical point. A scaling form for this 



function (valid for large N) may be written ji4 



p(r) oc r K exp(-Dr s ) (21) 



where S = (1 — v) 1 , D is a constant and k = 0.249 ± 0.011 in three dimensions ||44||. In 



figure P we plot the function p( r 2 ) r -°- 1245 against r 2 on a logarithmic scale, for the chain 
length N = 40 and N = 60. Fits to the data yields estimates v' = 0.51 and v = 0.518 
respectively. Of course it is hard to believe that the chains are swollen at T C (N) which is 
below 0, given the fact that for <p —>■ at T = 0, the chains are not swollen. It is thus 
possible that the slight deviation from u' — 1/2 is simply due to corrections to scaling. 



IV. DISCUSSION 

In summary we have performed a study of the liquid-vapor critical point of a polymer 
model for chain lengths up to N = 60 monomers. Owing to the low acceptance rate for chain 
transfers it was not possible to study either very long chains, or very large systems. Nev- 
ertheless we believe that the FSS-based technique we employed, of matching the measured 
scaling operator distribution functions to their fixed point universal forms is considerably 
more reliable than the practice of simply extrapolating a power law fit to coexistence curve 
data obtained well away from criticality, as has been the norm in previous simulation stud- 
ies p3| , p4] , p7| , p5]| . Indeed, while we reproduce the previous results with regard to the finding 
that T C (N) is well described by a Flory formula, our measured value for the exponent 
Xi = 0.37(2), is in much closer accord with experiment (X2 = 0.38) than previous simulation 
measurements. It is also interesting to note that for a similar range of monomeric units 
N, our coarse grained model seems to be much better at describing the asymptotic limit 



than chemically realistic models such as that employed in a the recent study of Alkanes p3 
which did not even yield a monotonically decreasing <j) c (N). 

With regard to the critical iV dependence of the chain span, our results suggest (albeit 
on the basis of a very limited number of data points) that the exponent v' > 0.5, which, 
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if correct, would imply that the chains are slightly swollen at criticality — at variance with 
the suggestion of Cherayil [19.46] and L'huiller [2(J that v' « 0.46, which is based on the 
assumption that the phase separation occurs when the chains just barely begin to overlap. 
We consider it possible, however, that our estimates for v' should be considered as effective 
exponents, which exceed the classical value v' = 1/2 only because of corrections to scaling. 
But in any case there is no evidence that the chains are somewhat collapsed at criticality. 

Finally we remark that there is evidently a need to study longer chain and larger system 
sizes in order both to validate the results thus far obtained and to confirm that the limiting 
scaling behaviour is being observed. In view of the low acceptance rates for chain transfers 
at large N, algorithmic improvements are clearly necessary before this can be achieved. In 
this regard, recent improved biased growth techniques for chain insertion at large iV and <p 
promise to be extremely helpful In future work we hope to report on their application 
to the present problem. 
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FIGURES 




FIG. 1. Schematic phase diagram of a polymer solution in the space of the temperature T 
and the volume fraction <j>. The ceoxistence curve separates a dilute solution of collapsed chains 
[at </>coex] from a semildilute solution of overlapping chains [at (0}ex\- These two branches of the 
coexistence curve merge at a critical point T C (N), (p c (N). For N — > oo the critical point merges with 
the G-point of a dilute polymer solution [T C (N — > oo) — > G, (f> c (N — > oo) — > 0] and the unmixing 
transition has a tricritical character. At T = G, the chain configurations are ideal Gaussian coils, 
while their structure at T C (N) is non-trivial. 
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FIG. 2. The CBMC acceptance rate as a function of the monomeric volume fraction for chain 
lengths N = 20, 40, 80. 



13 



15.0 



10.0 







5.0 



0.0 



N=64 




N=80 




N= 1 00 




N=150 


yy>^ — 






✓ / ' 

/ / 

/ / 






i 


i 



1.5 



1.8 2.0 

Temperature 



2.2 



2.5 



FIG. 3. The temperature dependence of Rg/N for a single chain with a variety of chain lengths. 
All simulations were performed at the temperature T = 2.0 and the temperature dependence ob- 
tained by histogram extrapolation. The common intersection point for large N yields the estimate 
6 = 2.020(3). 
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FIG. 4. (a) The distribution function of the monomeric volume fraction pl{4>) for the N=20 
system at a selection of temperatures along the line of liquid-vapor coexistence, (b) The measured 
liquid (squares) and vapour (circles) peak densities corresponding to the histogram extrapolation 
along the coexistence curve. Also shown is an attempted fit to the data (full line) and the estimate 
for the critical point (X), both obtained as described in the text. 
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FIG. 5. The ordering operator distribution pl(-M) of the polymer model at the assigned critical 
parameters. Also shown for comparison is the universal fixed point form pj^(x) obtained in a 
scpcrate study ]3{|. In accordance with convention all data has been scaled to unit norm and 
variance, by choice of the scale factor aj^. 
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FIG. 6. The measured estimates of the inverse critical temperature as a function of inverse 
chain length, together with the G temperature, representing the critical temperature in the limit 
of infinite N. Also shown is a fit of the form 1/T C = 1/0(1 + 0.402iY-°- 5 + 0.696N- 1 ). 
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FIG. 7. The measured estimates of the critical density as a function of chain length, together 
with the infinite N value <j) c = Also shown is a fit of the form <j) c = (1.1126 + LSiY - 369 )" 1 
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FIG. 8. The N dependence of the end-to-end distance squared. The fit is of the form 
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FIG. 9. The measured function p{R^,)R e °- 1245 . Fits shows that v > 0.5 implying that the 
chains are slightly swollen at criticality. 
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